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Abstract 

We study a non-perturbative formulation of J\f = 4 super Yang-Mills theory (SYM) 
on R x 5" 3 proposed in arXiv:0807.2352 . This formulation is based on the large- iV 
reduction, and the theory can be described as a particular large- iV limit of the plane 
wave matrix model (PWMM), which is obtained by dimensionally reducing the original 
theory over S 3 . In this paper, we perform some tests for this proposal. We construct an 
operator in the PWMM that corresponds to the Wilson loop in SYM in the continuum 
limit and calculate the vacuum expectation value of the operator for the case of the 
circular contour. We find that our result indeed agrees with the well-known result first 
obtained by Erickson, Semenoff and Zarembo. We also compute the beta function at 
the 1-loop level based on this formulation and see that it is indeed vanishing. 
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1 Introduction 



The AdS/CFT correspondence [IH3] has been intensively studied and fruitfully extended to 
various directions over a decade. However, a complete proof of this conjecture is still missing 
even for the most typical example, the correspondence between J\f = 4 super Yang-Mills 
theory (SYM) and type IIB superstring theory on AdS^, x S 5 . This is partially because the 
correspondence is a strong/weak one. Namely, the region on the string theory side in which 
the supergravity approximation or the classical string approximation is valid is mapped to 
the strongly coupled region in the planar limit on the gauge theory side. In order to study 
M = 4 SYM in the strongly coupled region, one needs to have a non-perturbative formulation 
such as the lattice gauge theory. In fact, there are considerable developments in the lattice 
theories of M = 4 SYM [4T4TT]. However, any lattice formulations of M = 4 SYM proposed 
so far seem to require 

In [12], a non-perturbative formulation of M = 4 SYM on R x S 3 in the planar ('t 
Hooft) limit was proposed by using the plane wave matrix model (PWMM) [13] (for earlier 
discussions, see [HE)])- Note that M = 4 SYM on R A at a conformal point has the 
PSU(2, 2|4) symmetry with thirty-two supercharges and is equivalent to M = 4 SYM on 
R x S 3 through the conformal mapping. The PWMM can be obtained by dimensionally 
reducing M = 4 SYM onRx S 3 over S 3 [THlfTTj . In the formulation, J\f = 4 SYM on R x S 3 
is retrieved by taking an appropriate large- N limit for the theory around a particular vacuum 
of the PWMM. Thus the formulation is viewed as an extension of the large- N reduction [18], 
which asserts that the planar limit of a gauge theory is described by its dimensionally 
reduced model. The formulation provides a matrix regularization of the planar M = 4 
SYM on R x S 3 , where the matrix size plays a role of the ultraviolet cutoff. Remarkably, 
the formulation preserves sixteen supercharges and the gauge symmetry. Note that this 
number of the preserved supercharges is optimal, because any regularization must break 
the conformal symmetry so that the number is inevitably reduced to less than or equal to 
sixteen from thirty-two. Since the formulation preserves so many supercharges and provides 
a massive theory, no fine-tuning should be needed in taking the continuum limit. This is 
advantageous to the lattice formulations of M = 4 SYM, while the latter may also be used 

1 It was shown recently that only a fine-tuning of one parameter is needed at least to one-loop order in a 
lattice formulation [11] . It is discussed that no fine-tuning is required in the recent proposals [9]ll0j. 



1 



to study the finite N M = 4 SYM. Thus the formulation gives a feasible way to simulate the 
planar A/_= 4 SYM by computer using the methods in USED] and study its strong conpling 
dynamicaj. 

To check the validity of the formulation, some one-loop perturbative calculations were 
performed in [T2]. It was shown that the tadpole vanishes and that the fermion one- loop 
self-energy agrees with the one in the continuum theory. In [271128] , the known result of 
the confinement- deconfinement transition in the planar M = 4 SYM on R x S 3 in the 
weak coupling limit at finite temperature was reproduced in the formulation. In [2H], a 
two-loop calculation in the high temperature limit was done, and the result was consistent 
with the continuum theory. In [30l|3T], based on the earlier work [32l|33], the same large- N 
reduction for S 3 was applied to Chern-Simons theory on S 3 to obtain a matrix regularization 
of the theory in the planar limit, and the exact results on Chern-Simons theory on S 3 
were reproduced. Extension of the formulation to other supersymmetric gauge theories was 
discussed in J3HES]- In [2S1E7J, large- N reduction on general group manifolds and coset 
spaces was found. Applying the case of SU{2) ~ S 3 to planar Af = 4 SYM on R x S 3 
yields another non-perturbative formulation of the theory in terms of the PWMM. Recently, 
a matrix model regularization of M = 4 SYM on S* 4 was proposed in |38j . 

In this paper, we make a further check by calculating two quantities: the vacuum expec- 
tation value (VEV) of a half-BPS Wilson loop and the beta function. 

Erickson, Semenoff and Zarembo gave a successful example of calculation in the strongly 
coupled regime of A^ = 4 SYM (see also [ID])- They considered (locally) supersymmetric 
Wilson loops in M = 4 SYM on R A . It is conjectured that the logarithms of the VEV of 
such a Wilson loop in M = 4 SYM corresponds to minus the area of the minimal surface of 
a string world sheet in the AdS$ x S 5 such that the boundary of the world sheet coincides 
with the loop [4T1I32]. In the planar limit, they evaluated the VEV of the circular Wilson 
loop, which is half-BPS, by summing up all ladder diagrams of all orders in the perturbative 
expansion. They indeed showed that in the strong coupling region the result agrees with the 
prediction on the gravity side [HHH]. Thus, they gave a nontrivial check of the AdS/CFT 
correspondence. Their result was reproduced through the localization method in [45] to 

2 Monte Carlo simulation of the correlation functions of the chiral primary operators in Af = 4 SYM and 
comparison of the results with the prediction from the gravity side will be reported in [21 j - Preliminary 
results including the one for the Wilson loop are seen in [22 j l23 j . For related work, see [24Tl26| . 
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show that it is indeed the exact result. In this paper, we reproduce their result in the above 
formulation of M = 4 SYM on R x S 3 . We first construct an operator in terms of the 
matrices in such a way that it coincides with the Wilson loop in M = 4 SYM on R x S 3 in 
the continuum limit. Then, we concentrate on the circular loop, and calculate its VEV by 
summing up all the planar ladder diagrams. We find that the VEV exactly agrees with the 
result in [39] . 

We also evaluate the 1-loop beta function and find that it indeed vanishes, which is 
consistent with restoration of the superconformal invariance. 

This paper is organized as follows. In section 2, we review the non-perturbative formula- 
tion of the planar M = 4 SYM on Rx S 3 proposed in [12] . In section 3, after reviewing Wilson 
loops in the AdS/CFT correspondence, we construct operators in the non-perturbative for- 
mulation which correspond to the Wilson loops in M = 4 SYM on Rx S 3 in the continuum 
limit. In section 4, we calculate the VEV of the operator in the half-BPS case and show that 
it agrees with the known exact result. In section 5, we calculate the 1-loop beta function 
and verify that it indeed vanishes. Section 6 is devoted to summary and discussion. In 
appendices, some details are gathered. 

2 Large- N reduction for M = 4 SYM on R x S 3 

In this section, we review the novel large- iV reduction for M = 4 SYM on R x S 3 proposed 
in [T2]. In section 2.1, as a warm-up, we review the large- N reduction on S 1 developed in [T2] , 
by taking the 3 theory on S 1 as an example. In section 2.2, to explain the mechanism of 
the large- N reduction on S 3 , we consider the large N reduction for the <j) 3 theory on S 3 . In 
section 2.3, we apply the mechanism to M = 4 SYM on R x S 3 . 

2.1 Large- N reduction on S 1 

We consider the 3 matrix quantum mechanics at finite temperature with the inverse tem- 
perature R. In other words, it is a one-dimensional matrix field theory on S 1 with the radius 
R. The action is 

*-?r"(K2)"+£*^)- 

where < x < 2nR, <ft is an iV x iV hermitian matrix, and £ is a dimensionless mass. 
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To diagonalize the quadratic part in (12. ip . we make the Fourier expansion 



(2.2) 



The action (12. ip is rewritten in terms of the Fourier modes as 



S 



Vs 
9 2 



2 ^ & 2 + £ 2 )^- + ^ E ^ 



+n2+n3,0 ( Pni ( Pn2 ( Png 



where V51 = 2ttR is the volume of S . 

Here, for instance, we consider the following correlation function 

— Tr(0(xi)0(x 2 )) 



(2.3) 



(2.4) 



By using the translation invariance of the theory and the Fourier expansion (12. 2p . we see 
that the above correlation function is equal to 



2^-R J d y (JjTrMx! - x 2 + y)<f>{y))^ = ^ (^ Tr ( 



g if (xi-x 2 ) 



(2.5) 



We calculate a diagram in Fig. 1 which appears in the perturbative expansion of (-iTr(0 n 0_ n )) : 



-Yy— 



1 



v s i ) " w+ej ^p+e(n+i) 2 +e 

This survives in the planar limit where 

N ->■ 00 with # 2 iV fixed . 



(2.6) 



(2.7) 



In this limit, all the planar diagrams such as Fig. 1 contribute to (12 ,4p . while other non-planar 
diagrams such as Fig. 2 do not. 

To obtain the reduced model of (12. ip . we introduce a constant diagonal matrix with the 
eigenvalues uniformly distributed: 



R V 2 ' 



2 ' ' 2 



where v and k are integers satisfying^ 



N = vk 



(2.8) 



(2.9) 



3 The tensor product of lfe is needed for the theory on S 1 to extract the planar contribution through the 
k —> 00 limit [T2]. It is not needed for the theory on R. 
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Figure 1: Planar diagram 
for the two-point correla- 
tion function 



Figure 2: Non-planar dia- 
gram for the two-point cor- 
relation function 



v turns out to play the role of the UV cutoff. The rule for obtaining the reduced model is 

(f)(x) e lPx <f)e- iPx , V s i v , (2.10) 

where is an N x N hermitian matrix independent of x. v is determined shortly such that 
the reduced model reproduces the original theory. The rule (I2.10p implies that 

^ie iPx [P,d>]e- iPx . (2.11) 

ax 

Then, applying the rule (I2.10j) to ( 12. ip yields the reduced model of ( 12.1 ft 

Sr = ^(_i[P, tf + |_ 02 + _L^ . (212) 

If is decomposed to a tensor product of a z/ x z/ matrix and a /c x k matrix following ( \2.8\\ . 
(I2.12p is expressed as 



S r = — tr 
T 



5> - ^ + f))0 (s 'V M) + ^ £ ^ W M ) , (2-13) 



where is a fc x matrix, and tr stands for the trace over k x k matrices. The range of 
the indices s, t, u is 



v + 1 v — 1 , 

<s,t,w< , (2.14) 



2 2 
and the indices s, t, u run integers for odd v and half-integers for even v . 

By applying the rule (I2.10p . we obtain the observable corresponding to (12.41) 

jjTi{e iPxi <f)e- iPxi e iPx2 <j)e- iPx2 )\ , (2.15) 



where (• • • ) r stands for the VEV in the reduced model. By using (I2.8p . this is rewritten as 

I^/ltr(0( s 'V ,s) )) e^ 1 " 3 ^ . (2.16) 



s,t s ' r 

We calculate again the diagram in Fig. 1 which appears in the perturbative expansion of 

(itr(0( s -*VC' s ))) r : 

V f^Yi? 6 ( 1 V 1 1 (2 17) 

) \{s-t) 2 +e) {t- u ) 2 +e{u-s) 2 +e ■ l ' ; 

Comparing (I2.6P and (12.171) . we put 

Vgx 

v = — 

V 

and take a limit 



(2.18) 



k -»• oo, v -> oo with g 2 N fixed . (2.19) 

Note that (I2.19P implies (12 .7p because N = vk. We make an identification n = s — t, 
I = t — u and — n — I = u — s , and find that (12.171) agrees with (12 .6p in the limit fl 2 . 1 9 [) . 
This correspondence holds for all the planar diagrams such as Fig. 1. As for the non-planar 
diagrams such as Fig. 2, such correspondence does not hold. However, because of k — > oo, 
only the planar diagrams contribute to ( 12. 161) . Furthermore, we compare ( 12.16P with ( 12. 5p . 
While the summation over s and t in (12.161) is redundant compared to the summation over n 
in (12. 5p . this redundancy is canceled by a factor 1/z/ in (I2.16p . Thus we conclude that (12.151) 
in the limit (12.191) agrees with (12.4p in the planar limit. 

In a similar manner, one can easily show that in the limit (I2.19p 

j^Tr(e iPxi (f)e- iPxi e iPx2 0e- iPx2 ■ ■ ■ e iPxr <f)e- iPxr )\ = /-^Tr(0(xi)0(:z 2 ) ■ • -<f)(x r )) 

(2.20) 

One can also show the following equality for the free energies of both theories in the limit 

F F 

(2.21) 



N 2 v N 2 V s i 

This gives an interpretation of v that it is the volume of space on which the reduced model is 
defined. Thus, in the limit ( I2.19p . the reduced model (" 12 . 1 2[) with ("12. 18[) retrieves the planar 
limit of the original theory (12. ip . 
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2.2 Large-iV reduction on 

To illustrate the large- N reduction on S 3 , we consider the <fi 3 matrix field theory on S 3 : 
Here the radius of the S 3 is 2//z, and its volume Vs3 is given by 



V, = ^ . (2.23) 

represents the volume element of unit S 3 . 4>{£lz) is an N x N hermitian matrix scalar 
field on S 3 . Ci is the Killing vector on unit S 3 . Note that S 3 is identified with the SU(2) 
group manifold. From this point of view, dQ 3 is the Haar measure which is left and right 
invariant, and Ci is the generator of the left translation on the group manifold. For an 
explicit form of dQ 3 and Ci, see (1A.7I) and flA.llj) . respectively. £ is a dimensionless mass as 
before. 

In what follows, the argument goes parallel to the one in the previous subsection. Indeed, 
S* 3 is an S^-bundle over S 2 , and this S 1 is a counterpart of the S 1 on which the <j) 3 theory is 
defined in the previous subsection. 

To diagonalize the quadratic part in (12.221) . we expand <fi in terms of the spherical har- 
monics on S* 3 defined in flA.121) : 

j 

0(fi 3 ) = ^ } j (f)j m mXjrnm{^3), (2.24) 
J m,m=~J 

where J take non-negative integers and half-integers. The properties of the spherical har- 
monics on S 3 are summarized in appendix A. m is viewed as the momentum along the S 1 
fiber p^2lfT^ rT5] . By using them, the action (12.221) is rewritten in terms of the modes as 



s v - 



9 2 



tL (-1)"( J( J + 1) + e 2 )Tr(0 Jmm 0j_ m _ m ) 

Jmfh 

" (2.25) 



C-J 1 m\rh\ J-imifhi J^m^m^^ijt' 'Jimifhi^p ^ J^ra^m^ 'J^m^fh-j, 

J\m\fh\ J2iri2fh2 J^m^fh^ 

where C Jimimi j 2 m 2 m 2 j 3 m 3 m 3 is defined in (IA.18|). 
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By identifying the S 3 with the SU(2) group manifold, we denote 0(^3) by <f)(g), where 
g G SU(2). We consider the following correlation function which is a counterpart of (12. 4ft : 

i-Tr(0(^)0(^))V (2.26) 



N 

The invariance of the theory under the right translation implies that (12.261) is equivalent to 

/ S {^<^i92 l 9)m))- (2-27) 
We substitute ( I2.24p into the above quantity: 

J 2^ lj^^{4'Jmm < Pj'm'm')) -Rj(fi , 25'l l )m,''mXjm'' fh{^z)Yj'm'm,' iS^z) 
JmmJ'm'm' ' 

= (-l) m H(MrVA(^%™W). (2-28) 

Jmfhm' 

where Rj(g) is the representation matrix of g in the spin J representation, and we have used 
(IA.12P and (IA.16p . We have obtained a counterpart of ( 12. 5p . We calculate a diagram in 
FigJH which appears in the perturbative expansion of (^Tr(0j mrfi 0j m '-m)) : 

/_-Qm+m' /g^NV ( 1 \ 2 (-l)™2-™2 (_l)m 3 -m 3 

J mm mJ2 mi Jz rri3 "13 

X Cj—m— m J2.m2.rh2 J r smo,m-fij— m'm J^—m^ — 7713 J?.mi—mi ■ (2.29) 

This survives in the planar limit (12. 7p . As before, in this limit, all the planar diagrams such 
as FigJH contribute to (12. 26j) . while other non-planar diagrams such as Figj5]do not. 

Next we construct the reduced model of (12.221) . For this it is convenient to introduce an 
A^-dimensional reducible representation of SU(2) in which the generators take the form 



(1/-1J/2 

e 

=-("-l)/2 



® 1*.) , (2.30) 



which represents multiple fuzzy spheres. Here are the spin j representation matrices of 
the SU(2) generators obeying 

[Lf,Lf)=ie ijk Lf. (2.31) 
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v is the number of different irreducible representations of SU(2), and the index s runs 
integers for odd v and half-integers for even v. The parameters j s and k s in (I2.30p satisfy 
the relation 

(u-l)/2 

E (2j s + 1)K = N. (2.32) 

s=-(j/-l)/2 



In particu 
parameters 



ar, in order to construct a reduced model of (I2.22p . we choose the following 



2j s + 1 = n + s, ^ = fc. (2.33) 
n is the centered value of the dimension of the representations, so that 

N = nvk. (2.34) 
The rule for obtaining the reduced model, which is a counterpart of (I2.10p . is 

<f>(g) ->■ G~ x <j)G, V s s -> v, (2.35) 

where G is the representation matrix of g in the reducible representation specified by (I2.33p . 
and <j) is an N x iV hermitian matrix independent of f2 3 . v is determined shortly such that 
the reduced model reproduces the original theory. By using (1A.9[) . it is easy to see that 

d(j) -»• G~ l [Li, 4>]G. (2.36) 

Then applying the rule to (I2.22p yields the reduced model of f )2.22p 

Sr = ^Tr {-^[UAf + + |0 3 ) ■ (2.37) 

To see how this reduced model retrieves the original theory, we expand the (s,t) block 
of 4> we denote by (j)( s,t > in terms of the fuzzy spherical harmonics defined in (1B.3[) : 

<t> M = E E (2-38) 

J=\js-jt\ m=-J 



4 The correspondence between the notation here and that in [12] is (^here = (^o)there' Where = 
+ ^)there' (^)here = (^ s )there' (^)here = (^)there- 
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where is a k x k matrix. The properties of the fuzzy spherical harmonics are summarized 
in appendix B. We will see shortly that j s — j t = (s — 1)/2 is identified with rh as in the case 
of the 3 theory on S 1 . Using the modes, f!2.37j) is expressed as 



S ,■ 



vn 



t E E(- 1 ) m " 0s ~ it) ( J ( J + !) + ^M&^Stt) 



s,t Jm 



q E E 



(jsjt) J2m 2 (j t ju) Jamaijujs^W Jxmx^ Jim-i? J 3 m3 



i (s,t) i ; («,s) 



(2.39) 



where Cj imi (j a j t \ j 2 m 2 (jd«) J3m 3 (j u j,) is defined in (IB.8|) . and tr stands for the trace over /c x fc 
matrices. 

Applying the rule (12.351) . we obtain the observable in the reduced model corresponding 
to fl2T26|) : 

1 



N 



(2.40) 



where (• • -) r stands for the VEV in the reduced model. By using (I2.38p . flB.4|) and flB.6j) 
we calculate (I2.40D : 

s,t JmJ'm' ' ' ' 

;/ > > (-1)"' -••"-'-7,',( /A >.vr t )-,/»,( 1 



m"m 11 \Xjm" (jsjt)^J'm' (jtjs) > 



s,t Jmm' ' ' r 

We calculate again the diagram in FigJH which appears in the perturbative expansion of 



E E 



-Qm+m' / g 2 j\T \ 2 / | \ 2 ('_1'|m2-(i s -Ju) C_l^m3-Cj'«-ii) 



ii Jmm' J 2 m2Jzm 3 



nVvy v^ + i) + £V ^(J 2 + i) + e 2 ^(^3 + i) + e 2 



x Cj-m(jtjs) J 2 m 2 (j s j u ) J3m'i(jujt)Cj-m'(j s jt) Js-m 3 (J t j u ) J 2 m 2 (j u j a )- 



Comparison of (I2.29p with ( I2.42p leads us to put 

V S 3 



9 ' 



(2.42) 



(2.43) 



and take a limit 



k — > oo, f — ?• oo, n/v — > oo with g 2 iV fixed. 



(2.44) 
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Note again that (I2.44p implies (12. 7p . Then, in the limit (I2.44p . we can make an identification 
fh = j s — j t , m 2 = is — ju, ™3 = ju — jt as anticipated and show by using (IB.9P that 



Cj-m(jtjs) J2m2(j s ju) J'im-i{j u jt) ~ ^ Cj-m-m J2m2VH2 Jsm^ms 
Cj-m'(jsjt) Jz-m^(Jtju) J2m,2(jujs) ^ ^ J—m'rh Jz-ms-mz J2rn2-fri2 • (2-45) 

We therefore find that (I2.42p agrees with (12.291) . As before, this correspondence holds for all 
the planar diagrams such as FigJH while it does not hold for the non-planar diagrams such as 
Figj5] However, because of k — > oo, only the planar diagrams contribute to (I2.4ip . Note also 
that the UV/IR mixing on fuzzy spheres does not exist in the planar contribution, so non- 
commutativity does not remain in the n — > oo limit. Furthermore, while the summation over 
s and t in (I2.4ip is redundant compared to the summation over fh in (12.281) . this redundancy 
is canceled by a factor \jv in (12.411) . Thus we conclude that (12.401) in the limit (I2.44p agrees 
with (I2.26P in the planar limit. 

As in the case of the 3 theory on S 1 , one can easily show that counterparts of (12.201) 
and fl2T2T]) hold in the limit (12^) : 

^Tr(GrVGiG 2 VG 2 ---G T TVa)^ = ^Tr(0(^)^ 2 ) • • • 0^))^. (2.46) 

and 

= F (2 47) 

N 2 v mVss' 1 ' J 

Thus, in the limit (I2.44p . the reduced model (I2.37P with (I2.43P retrieves the planar limit of 



the original theory (I2.22p l 5 



What we have done is understood as follows. 5* 3 is viewed as an S^-bundle over S* 2 . 
The momenta along the S 1 is given by fh in (I2.24p . By making the Kaluza-Klein (KK) 
reduction along the S , one obtains the KK modes on S 2 with the KK momenta given by 
fh. The KK mode with the KK momentum fh behaves as in a situation that the monopole 
with the monopole charge fh is located at the origin of the S 2 . Reflecting the angular 
momentum possessed by the monopole, the angular momentum J of the KK mode on S* 2 is 
restricted to \fh\ < J. The identification fh = j s — j t = (s — t)/2 matches the range of J in 



5 Another large- N reduction on S* 3 was developed in [36], where the planar limit of the original theory is 
retrieved by setting n = (v + l)/2 and k s = (2j s + l)k in (|2.30[) instead of (|2.33|) and taking the limit in 
which v — > oo and k — > oo with v — Vs 3 k/N. 
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(1 2 . 3 8 [) in the n — > oo limit. Actually, the monopole spherical harmonics for the monopole 
charge q is regularized by the fuzzy spherical harmonics which is an (I + q) x / rectangular 
matrix [12 | [T4 ^ [T5] . The relation (I2.45P reflects this fact. We obtain the base space S 2 through 
the continuum limit of the fuzzy sphere and the S 1 fiber through the mechanism of the large- 
N reduction on S 1 explained in the previous subsection, n plays the role of the UV cutoff 
on S 2 , while v plays the role of the UV cutoff on S 1 . 

2.3 M = 4 SYM on R x S 3 from PWMM 

We first see that J\f = 4 SYM on R at a conformal point, where the VEVs of all the scalar 
fields vanish, is equivalent to M = 4 SYM on R x S 3 through the conformal map. In this 
subsection, for simplicity, we suppress the terms including fermion fields. Then, the action 
of M = 4 U(N) SYM on R 4 takes the form 

Sym = ^- 2 J d'xTr (f% + ^(D^ m ) 2 - \[<t>m, 0n] 2 ) , (2.48) 

where _D M = <9 M + ]. There are six scalars 4> m , and for later convenience we make m 

and n run from 4 to 9. This theory possesses the PSU(2,2\4) superconformal symmetry 
with 32 supercharges. Applying the Weyl transformation defined by 



Aft — >■ , 

4>m -> e 2 <p m , 



V -> ^ = e'(%„ . (2.49) 
to (12.481) yields Af = 4 SYM on a curved space endowed with a metric g^: 



(2.50) 

where R is the Ricci scalar constructed from g^ u . i? 4 is transformed to R x S 3 by the Weyl 
transformation. In fact, if one starts with the metric of R 4 in the polar coordinate and 
rewrites it as 

ds 2 Ri = dr 2 + r 2 dQ 2 3 = e" T (dr 2 + (^j dQ 2 ] , (2.51) 
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one obtains the metric of R x S 3 up to the conformal factor e MT . Namely, p = —fir in this 
case. Here we have changed the coordinate from r to r = ^ log (^) where 2//x is the radius 
of the resultant S 3 as we have adopted so far. Thus M = 4 SYM on R A at the conformal 
point is equivalent to Af = 4 SYM on i? x S 3 . 

Next, using the equations in appendix A, we rewrite the action (12.501) for R x S 3 in terms 
of the Killing vector. We expand the gauge field on S 3 as 

A = AV, (2.52) 



where e l (i = 1,2,3) are defined in (1A.3I) . and put m = X m . By using (1A.5I) . we rewrite 
(I2.50p in terms of Aj and A m : 

2 



S- 



= ^ J dr^Tr (~{D T Xi - ifxCiA T ) 2 + ~ (pXt + te ljk (fxC 3 X k + ~[A i; A fc ])^ 

+ -(-D T A m ) 2 — -(/i£jA m + [Aj, A m ]) 2 + ^r-^m — ^[A m , A n ] 

(2.53) 

where D T = d T + i[A T , ]. Applying the rule ( 12. 35ft to the above action, we obtain a reduced 
model of M = 4 SYM on R x S 3 : 



q - V 

JYM,r — — 



J drTr (^(D T Xi-in[Li,A T }) 2 + ^ (pXi + ie ijk {^[L h X k \ + ^[X 3 -, X^ 

+ l -{D T X m f - ^[L h X m ] + [Aj, A m ]) 2 + ^X 2 m - l[X m , X n f 

(2.54) 

Note that the matrices still depend on r. 
Here we can absorb /iLj into Xf 

fxLi + Xi^Xi. (2.55) 

In other words, we can regard /iLj as a background of Aj. Note that one can not perform 
such an absorption in the case of the 3 theory in the previous subsection. The resulting 
action is nothing but the plane wave matrix model 



Spw — — — / drTr ( -(D T X M ) 2 — -[X M ,X N ] 2 + ^—X 2 + ^—X^ + ipe i3k Aj Xj X k J 
9pw J \ z 4 1 ° J 



(2.56) 
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where M runs from 1 to 9 and we put 



9pw 



(2.57) 



v 



This action is obtained by dimensionally reducing M = 4 SYM on R x S 3 over S 3 [T6"|[l~7] . 

The PWMM possesses the SU(2\A) supersymmetry, which includes 16 supercharges. The 
PWMM possesses many discrete vacua given by Xj = fiL iy where Li is given in ( I2.30h . and 
they are degenerate. There are no classical moduli, and these vacua are classically stable. 
Furthermore, all of these vacua preserve the S77(2|4) supersymmetry, and the theories around 
these vacua are massive. Hence they are also quantum mechanically stable at least at 
perturbative level [16] . We pick up a particular one specified by f !2.33|) from these vacua and 
take the limit (I2.44p to obtain the reduced model of M = 4 SYM on R x S 3 . In this limit, 
this vacuum is stable even at non-perturbative level, because tunneling to other vacua by an 
instanton effect is suppressed by the k — > oo limit. The global gauge invariance (depending 
only on r) of the PWMM is translated to the local gauge invariance of M = 4 SYM. The 
SU(2\4) symmetry is expected to enhance to the PSU(2,2\4) symmetry in the continuum 
limit. Thus the reduced model (I2.54p retrieves the planar limit of TV = 4 SYM on R x S 3 as 
the reduced model (I2.37P does the planar limit of (12.221) . 

Note that obtaining a reduced model of M = 4 SYM on R 4 is quite non-trivial. For 
instance, dimensionally reducing it to R yields the DO-brane effective theory or the Matrix 
theory [47] . which is obtained by putting jj, — in (I2.56p . This theory possesses infinitely 
many continuous vacua, which are characterized by [Xm,X^\ = 0. All of these vacua 
preserve sixteen supercharges as well. However, there are one-dimensional massless fields 
around these vacua at least perturbatively, so whether each vacuum is stable at quantum 
level heavily depends on dynamics. Thus one may not naively expect that a theory around 
a vacuum where Xj are simultaneously diagonal and their eigenvalues distribute uniformly 
in R 3 retrieves the planar limit of M = 4 SYM on R 4 . 

We summarize the prescription of the large- N reduction for M = 4 SYM on R x S 3 . We 
expand the PWMM (I2.56P around the vacuum Xj = /iLj with (I2.30p and (I2.33p . and take 
the limit 



k — )• oo, v — )• oo, n/u oo with 



9p\yk _ 9 



2 N 



fixed. 



(2.58) 



n 
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Then, the planar limit of M = 4 SYM on R x S 3 with the 't Hooft coupling g 2 N is retrieved. 
The correspondence is explicitly given in (I2.46P and ( I2.47p . Thus the planar limit of M = 4 
SYM on R x S 3 is regularized by the PWMM. This regularization preserves the SU{2\4) 
symmetry including sixteen supercharges and the gauge symmetry. This number of the 
preserved supercharges is optimal. In the next section, we will see the correspondence 
between Wilson loops. 

3 Wilson loops 

In this section, we examine the correspondence between Wilson loops in M = 4 SYM and 
its reduced model we reviewed in the previous section. In section 3.1, we review the Wilson 
loops considered in the gauge theory side in the context of the AdS/CFT correspondence. 
In section 3.2, we construct the corresponding Wilson loop in the reduced model. 

3.1 Wilson loop in the AdS/CFT correspondence 

The Wilson loop considered in the AdS/CFT correspondence takes the form 

W(C) = ^TrPexp U ds{x»{s)A^x{s)) + *|x( S )|0 m ( S )0 m (x( S ))}^) , (3.1) 

where the function x^(s) : [0, 1] — > C specifies the contour C, and m (s) satisfies 5 mn Q m Q n = 
1. The Wilson loop is invariant under the Weyl transformation (I2.49p . namely takes the same 
form both in R 4 and RxS 3 . The contour C on i? 4 is mapped to the corresponding contour C 
on Rx S 3 given by the coordinate transformation used in f)2.5ip . and the following equality 
holds: 

(W(C)) Ri = (W(C')) RxS s. (3.2) 

We next consider the supersymmetric property of the Wilson loop on R x S 3 for later 
convenience. Its infinitesimal variation under the supersymmetry transformation is propor- 
tional to (x M r M + i\x\Q m T m )e(x) where /i = r,9,(p,i/) and e(x) is a ten-dimensional Killing 
spinor on R x S* 3 satisfying the Killing spinor equation |49|l50]. 

V a e{x) = ±^T a T l T 2 T 3 e{x). (3.3) 
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Here a = r, 1, 2, 3 is the local Lorentz index. This equation is solved by 

/ \ Er Tl2S T i \ _iir Tl23 -r _Ir 12 n _Ir 31 fl _Ir 12 w> /r, ,\ 

e l (x) = e4 i T r] ll e 2 (x) = e * l T e ^ v e ^ °e ^ (3.4) 

where e\ and e 2 are the solutions to the upper and the lower signs of (I3.3p . respectively, and 
rji and r)2 are constant spinors. If the Wilson loop operator is invariant under some of the 
supersymmetries, there exist some non-zero components of ^,2 such that the following two 
equations hold, 

(x»rp + i\x\Q m r m )ei rr ™ r r )1 = 0, 

(AT„ + z|i|e m r m )e^ rTl2 ^e-5 rl2 ^-^ r3l9 e^ rl2 ^ 2 = 0. (3.5) 
Introducing local projection operators as 



P 1± = \e-^{l ± ^Q-T, m )e^ 
2 \x\ 



P 2± = \ e ^ e ¥^ e h^ e ^r {l ± !^e m r / , m )e-^ Tl23 ^-^ 12 ^-^ r31e e^ r12 ^, (3.6) 
one can rewrite ( 13. 5 p simply as 

Pi+?7i = 0, P 2+ ^ 2 = 0. (3.7) 

The number of independent non-zero components of 771 and 772 which satisfy ( 13. 7p (or equiva- 
lently ( 13. 5p ) is just the number of supersymmetries preserved by the insertion of the operator. 

As an example, let us consider the case where x^(s) and 0(s) are constant in s and the 
contour does not extend to the r direction. In this case, the projection operator Pi is a 
constant operator on the contour, so that the half components of rji which are projected by 
P\ + into can be non-zero constants. Therefore, the operator preserves at least 8 supersym- 
metries (1/4 BPS). We can also consider a special case of this example in which the path is 
given by a great circle on S 3 at a fixed value of r, 

C : (r(s), 6{s), (fi(s),i>(s)) = fa, 0, 0, 4vr S ), (3.8) 

where To is a constant. In addition to the above mentioned 8 supersymmetries, another 8 
supercharges are preserved in this case because P 2 is also a constant operator on the contour 
( 13. 8 p as well as Pi. Hence, the Wilson loop on the great circle ( 13 ,8p is a half-BPS operator. 
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The contour (13. 8p is mapped to the circular loop, whose center is located at the origin on 
R 4 , by the Weyl transformation (12.491) . For the circular contour with the unit radius, a; M (s) 
can be parametrized as 



{x^(s)} = (cos(27rs),sin(27rs),0,0). 



(3.9) 



Because of (13.21) . the VEV of the Wilson loop with the contour (I3.8P on R x S 3 coincides 
with the VEV of the Wilson loop with the contour (13. 9p on R 4 . 

The VEV of the BPS circular Wilson loop on R 4 can be computed including all orders 
in the perturbative expansion in the planar limit (39JSQ1 - Let us review the result of this 
computation. The VEV of (13. ip for the path (13.91) turns out to be given by the following 
VEV in the Gaussian matrix model, 

(^(circle)) = ^-Trexp(M)^ = ^ J VM^-Tr exp(M) exp ^-^TrM 2 ^ = ^hiV^X), 

(3.10) 

where A is the 't Hooft couplings and I\ is a Bessel function. The strong coupling expansion 
of this result reproduces the correct coupling behavior observed in the gravity side. 



3.2 Wilson loop in the reduced model 



By using the relation (I2.46p . we can construct an operator in the reduced model (12.541) 
corresponding to the Wilson loop (J37TJ) : 



W{C) = -Tr 



Pexpj* J ds (x t (s)G- 1 (s)A t (t(s))G(s) 

+ x*(s)e} l (xt(s))G- l (s)X l (T(s))G(s) 
+i\x(s)\Q m (s)G- 1 (s)X m (T(s))G(s))}] , 



(3.11) 



where fi = 6, ip, ip and e % is the right invariant 1-form on S* 3 defined in (IA.3[) . G(s) is given 
by 



G(s) = Pexp 



% \ ds'x^is'^ix^is'^fiLi 
o 



(3.12) 



3 Our convention of the 't Hooft coupling is slightly different from that in [39]. The relation is given by 
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where Li is given by (I2.30P and ( 12 .33j) . Note that G(l) = 1. By referring to an argument 
in [IB], we can show that (13. lip is equal to [32] 

Pexp ji J ds (x t {s)A t {t{s)) + x' 1 {s)e i fl {x p {s)){fiL i + Xi(r(s))) 
+i\x(s)\e m (s)X m (r(s)))}]. (3.13) 



W{C) = ±Tr 



In fact, (13. lip can be viewed as 1/N times the trace of the time evolution kernel for a 
time-dependent Hamiltonian 

H(s) = - {^( S )el(^( S ))/iL* + x t (s)A t (t(s)) + x fi (s)e i p (x fi (s))X i (T(s)) 

+z\x(s)\Q m (s)X m (r(s))}. (3.14) 

If one regards the first term as a free part and the others as an interaction part and switches 
the picture to the interaction one, one obtains (13. lip . Note that applying (I2.55P to (13.131) 
results in the dimensional reduction of (13. ip from R x S 3 to R. In other words, we expand 
Xi around /zLj in the dimensionally reduced Wilson loop. Such a correspondence for the 
Wilson loops is seen commonly in the large- N reduction. 
From (12.461) . we expect 

(W(C)) r = (W(C)). (3.15) 

In the next section, we verify this relation in the case of the half-BPS Wilson loop with the 
contour C given in (13.81) . 

4 VEV of the half-BPS Wilson loop from PWMM 

In this section, we calculate VEV of the operator (13.131) with the contour (13.81) . 

For the calculation of the circular BPS Wilson loop in SYM on R 4 , it is convenient to take 
the Feynman gauge where the computation is simplified such that the non-ladder diagrams 
cancel out [39]. One can expect that this simplification also occurs in the PWMM if one 
takes a particular gauge which corresponds to the Feynman gauge in the continuum limit. 
This leads us to take this particular gauge in the PWMM for the following computation. 

The corresponding gauge fixing term on R x S* 3 can be obtained by applying the Weyl 
transformation. Although the original action is invariant under the Weyl transformation 
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(I2.49p . the gauge fixing term is not invariant in general. In fact, if one takes the Feynman 
gauge on R 4 , the gauge fixing term is transformed as 

5^/5* (<«vf Af) 2 = {wVFtG** )} 2 



2V« 

e 



2/tT , . 2 



= (< xS aVf + 53 ) 2 • 

Hence, the corresponding gauge fixing term on Rx S 3 should be given by the last line of ( 14. ip . 
We then apply the reduction rule (I2.35j) to the above gauge fixing term plus appropriate ghost 
terms to obtain those in the PWMM, 

S gf+gh = — — / drTil -(d T A T + ifi[Li, X t ] + fiA T ) 2 
9pw J ^ 

— ic{d T D T + fiD T )c + fic[Li, ifi[Li, c] + i[c, j, (4.2) 

where c and c are the ghost fields. 

We add ( 14. 2 ft to the original action of the PWMM expanded around the background 
( 12.301) with the parameters given by ( 12.331) and construct the Feynman rule of this theory. 
We make a mode expansion of the matrices as in (1C.9j) and read off the propagators for each 
momentum mode as 

( x JmJ,(P)ii4w(p')«) 

^(-l^-^-^+Hjj^^Sp^S^SiiS^wS^ + P')^f (p + 0) 
9j^ { _ 1) m-Us-, t )+i 6jJl6m ^S^SuSj^Sip + i/) tf$$$$$ m (p = P ' = 0) 

2 

(B^WijB^Wu) = 9 -^{-l) m -^-^5jj l 5 m . m ,5 st ,5 ts ,5u5 jk 2^(p+p')- 



n ■ ■ (p 2 +M 2 J 2 )(p 2 +M 2 (J + l) 2 ) ' 



(c i jt\phc i J Z\p%) = &&(-rr-V'-M6jj,6 n „ m ,6 st ,6 ts ,6u6 jk 2ir6(p + p') _^ +i ^ J(J+1) , 

^ABJmiPh&j'rn^kl) = — ^ABA'B'(-l)^ W '- A) ^J'5 ro - m '^<W*«*i*27r*(p + j/)-^, 
(Piij^A' 'j'L'k'(p')m) = ^^SjJ'Smm'SKK'diiSss^tt'Siidj^irdip - p') P+ ""J , (4.3) 

where Uj, Uj and Co>j are defined in (1C.11I) . We can also read off the interaction vertices in 
this mode expansion. In particular the vertices used in this section take the same form as 



in AC. 131) although a different gauge is taken in appendix O 

19 



(T-a) (T-b) 




(T -d) (T- e) (T - /) 



Figure 3: Tadpole diagrams. The curly, wavy, dotted, solid and dashed lines represent the 
propagators of X iy A t , the ghost, $ab and ^ A respectively. 

In general, the tadpole is expected to vanish as in the continuum theory. Let us check 
it at the 1-loop level in the present gauge. The only possibly nonzero contribution is the 
truncated 1-point function for ^j^l(p)y, where i,j run from 1 to k and p is the momentum 
along the r direction. The 1-point function takes the form, 

27rd(p)d st d ij d p ^8 J0 8 m0 5 ij T^. (4.4) 

There are six diagrams for the 1-loop correction to (I4.4p as shown in Fig. |3j Note that all 
these diagrams are planar ones. The diagrams (T — a) and (T — b) completely cancel each 
other. The diagram (T — /) vanishes under the integration of p because it is an odd function 
of p. Below we list the value of Y^) for each of the remaining diagrams, 

(T - c) = h(-i) R+ ^WR(R + WR+i){l I (4-5) 

V n t,R 

(T-d) = 12^ Yl u-i) R+i ' +i wm+WR+i){l I (4.6) 

v t,R 

(T - e) = -16^| Y, h(-l) R+ ^VRjR + WR + ^{l f s f}, (4-7) 



t,R 

where the last factor in the above equations is the 6-j symbol. We therefore find that the 
1-loop contribution to the tadpole in the present gauge is indeed vanishing without taking 
the continuum limit. 

Let us calculate the VEV of the Wilson loop (13. lip with the contour given by a great 
circle (13. 8ft . Thanks to the SO{Q)r rotational invariance, we can set O m = 5 m4 without loss 
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of generality. Then, the operator (13. lip for the great circle can be written in a relatively 
simple form by substituting ()A.4j) for the circular contour, 



W = -^-Tr ( e 4m ^ +H A , 
vnk 

where we have defined a dimensionless complex matrix H as 

H=-(X 3 + tX 4 ). 
H is expanded in terms of the fuzzy spherical harmonics as 

■fii+lm „.(«>*) | y^Rm (s,t) _ ■ r iR-lm rr ,{s,t) 



RmlO X Rml 



RmlO^RmO ^ RmlQ X R- J ;„ - J. 



R,rn 
R,m 



Rm ^ 1 Rm > 



where we have used (IB. lOj) . (1C2[) and (1C.9|) . Here the second term 



(s,t) , 
RmO< 



(4.8) 



(4.9) 



(4.10) 



and the third 



term (x^l lm _ x ) in the right-hand side of the first line are summed over R>\ and R>1, 
respectively, while the first term (x^^) and the last term 



J AB,Rm, 



are summed over i? > 0. 



Using (14. 3 p and the explicit expression of the Clebsch-Gordan coefficients 

'( y R+l + m)( y R + l- m y 1/2 



/^iR+lm 
^RmlO 



/~iRm 

°i?ml0 [ n, r > i .11/2 



(2i2 + l)(i2 + l) 
m 



/~tR—lm 
^RmlO 



[R(R+1)Y/ 2 
(R + m)(R-m)] 1/2 



(4.11) 



i?(2i? + 1) 

we can calculate the propagator for H. In particular, the equal-time propagator is simply 
given by 



^Rm ( T )ij^R'm'' ( T )kl) ~ &R0 $R' ^m0 $m' O^st^s't'^st'^il &j 



r(s',t') 



9pw 



(4.12) 



2/rn 

Namely, only the zero mode contributes to the propagator of H in the present gauge. We 
expand the VEV of the Wilson loop as 

(W) 



jLv^l ( ^ Tr / e 47r»(L 3 +aJJ)A 

nk ^ q\ V dai 1 J „ n 

— \ / da i / d«2 ■ " • / da q+ \5{\ — a\ — a^ 
n k „_ n Jo Jo Jo 
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Figure 4: Planar ladder diagram Fi S ure 5: Planar non-ladder di- 

for Wilson loop agram for Wilson loop 



x (4ni) q Tr(e 4lTiaiL3 He 47ria2L3 H ■ ■ ■ e ^^Ee^ ia ^ lLz ) , (4.13) 

where we have used the generalized Feynman formula proven in appendix [D] to show the 
second equality. Then we obtain 

y 00 r-l r-l pi 

(W) = — - y. /, / doti I dot2 ■ ■ ■ da q+ i5(l — a\ — CX2 — ■ ■ ■ — a q 
unk q=0 {si} J ° J o Jo 

x ^nia^! e 4iria 2 r2 _ . . e ^-Kia q T q ^Kia q+Y r x ^| t ]\n -ji+r 2 -j 2 -\ \~r q -j q 

Y fiRimi (~iR 2 m 2 (~iR q m q ry / tt(s 1 ,s 2 ) tt(s2,s 3 ) Tr( s q, s i)\ (ATA') 

A Jifij2-f2 U j2f2j3-r3 " ' ^j g r q ji-n ±L \ II Rim 1 11 R 2 m 2 " ' 11 R q m q h K^-^) 

where C^!p r , is the Clebsch-Gordan coefficient. 

In the following, we assume that only the planar ladder diagrams (i.e. planar diagrams 
without vertices, see Fig. H]) contribute to the VEV as in the case of the continuum theory. 
This assumption is reasonable since we are working in the gauge corresponding to the Feyn- 
man gauge on R 4 . In the Feynman gauge, the ladder approximation gives an exact result in 
the case of the circular Wilson loop in SYM on R 4 . Hence, also for the PWMM, we believe 
that the non-ladder diagrams do not contribute to the VEV in the continuum limit in the 
present gauge. 

Within this assumption, only the modes with R = m = contribute to (W) because the 
tree-level propagator (I4.12p is vanishing unless R — m — 0. Furthermore, only the terms 
with q even contribute. Thus, we obtain, 
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where we have used C°^_ r = , 2 . +1 . Taking the limit (12.581) . we obtainL 

oo 

= 1 y XW-J-Trexp(M)exp ^-^TrM 2 ^ . (4.16) 

On the right hand side of the first line of (I4.16p . we have put (s, s) to (0,0) because the 
propagator does not depend on s. In the second line, we have put A = s ™ and identified 
47riifoo °^ with a k x k hermitian matrix M. Although H^' °' ) is a complex matrix and time- 
dependent, this identification is possible within the ladder approximation where only the 
equal-time propagator in the free theory (I4.12p is needed for the computation of the Wilson 
loop. Since the equal-time propagator of AniHQQ' ^ takes the same value as that of M if 
we take the weight of the integral of M as in (I4.16p . the computational rules are identical 
between the two expressions. The result (I4.16P agrees with (I3.10p . 

Here, we again emphasize that in a general gauge one needs to add also the planar non- 
ladder diagrams (See Fig. [5]) to obtain the correct result (13.101) in the continuum theory. In 
fact, we checked that in a different gauge which will be introduced in the next section the 
sum over only the planar ladder diagrams does not coincide with (13.101) . This is consistent 
with the fact that the value of each Feynman diagram depends on a gauge choice and the 
non-ladder diagrams cancel out only in the Feynman gauge on R 4 . 

5 Beta function 

In this section, we compute the beta function at the 1-loop level in the PWMM expanded 
around the background specified by (I2.33P and show that it is indeed vanishing. In this 
section, we take the gauge shown in appendix [C] for convenience which is different from what 
we used in the previous section. We consider only the planar diagrams in the following. 

We first compute the wave function renormalization of the SO(6) scalar field. The self- 
energy is given by the truncated two-point function ((t^ABjm^ij^cDJ'm'iP')^) which takes 
the form, 

2*5(p - J n8«8, t \ejw D W JJ ,8^*;>%). (5.1) 



Here, the limit v — > oo is not used explicitly. However this limit should be needed for the non-ladder 
diagrams (See Fig. [5]) to cancel out. 
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There are six diagrams contributing to the 1-loop correction to the self-energy as shown in 
Fig. O For example, one can compute (5 — 2) in Fig. [B]as follows. In terms of the Feynman 
rules shown in appendix [Cj the contribution of (S — 2) to Hj'^(p) in (15. ip can be written as, 

dq iq + KiwJ i(p~q) + K 2 UJ% ~j 2T 



- 2 *E E 

JimiKi J 2 m 2 K2 



2m 2 K 2 (jujt) 



Jl-miKi(j„j s ) 



27T g 2 — (yfl 3 ) 2 (p — g) 2 — (u)j 2 1p) 2 J i- m i K i(>>) Jm (^jt) J2m2K2(jujt) Jm(jijs)' 

(5.2) 



where the ranges of the variables in ^2 Jm . K . are given by «j = ±1, m 8 = — J;, — Jj + 1, • ■ ■ , Jj 
(z = 1, 2), Jx - «! = |j s - j n |, |j s - jj + 1, • • • , j fl + j u and J 2 - k 2 = \j t - jj, \j t - j u \ + 
I? ' ' ' ijt + Jw By performing the integration, substituting the explicit form of J 7 shown in 
(IB. 171) and finally taking the summation over m 1 ,m2, K\ and K 2 , one can obtain, 

js+ju jt+ju 

2>2ixnk{-l) m ~ {j °- ]t) ^ Yl (2J+l)(2 J Ri + l)(2i? 2 + l) 



Rl = \js-ju\ R2 = \jt-ju\ 



{R 1 + l){R 2 + l){R l + R 2 + l 



P 2 + /X 2 (.Ri + ^ + f 



3^2 



i?2 + 2 -^2 | 



J 



J 



i?l R\ 



1 1 

2 2 
1 1 



RiRajRi + Rz + \) j p p 



R 2 -Rl «^ 

Js Ju 



(5.3) 



where i?j are defined as Ri = Ji — Ki for = ±1. Then, if one sets J = and s = t — 0, 
it is easy to compute the divergent part of (15.31) . In this case, the explicit formula for the 
6j-symbol and 9j-symbol can be found in [52] • Substituting those explicit forms, one finds 
that the divergent part is given as 

(v-l)/2 n-l+u/2 

(5.4) 



^ u=-(i/-l)/2 R=\u/2\ ^ 

Similarly, the divergent parts in the other diagrams in Fig. [6] are evaluated as 

(f-l)/2 n-l+u/2 

( s -d4 E E <- 5 >. 



w=-(i/-l)/2 R=\u/2\ 
(v-l)/2 n-l+u/2 



< s - 3 > = ^ E E 



u=-(v-l)/2 R=\u/2\ 
(f-l)/2 n-l+u/2 



-3 + 



E E -sMH 

r «=-(i/-l)/2 fl=[u/2| k \t" / J 



\u/2\ 
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(5-4) 



(5-5) 



(5-6) 



Figure 6: The 1-loop corrections to the self-energy of scalar field. The solid, dashed, curly 
and wavy lines represent the propagators of the 50(6) scalar, fermion, 50(3) scalar and the 
gauge field, respectively. 

(f-l)/2 ti-I+u/2 - 
rv \ "V v "\ 



n u=-(i/-l)/2 R=\u/2\ 
k ^ "-^" 2 fl 1 C:,,r 



^ u=-(i,-l)/2 fl=[u/2| k V ^ 7 -* 

Here, it is important that the quadratic divergence cancels when one takes the sum over all 
the diagrams in Fig. [6j Note also that there is no divergence even if one takes the large-n 
limit. This limit corresponds to the commutative limit of fuzzy sphere in which the theory 
on R x 5 2 is realized. The fact that there is no divergence due to n is consistent with the 
super renormalizability of (2+1) dimensional gauge theory. 

From (15.4)) and ( 15.5)) . we find that the wave function renormalization for the 50(6) scalars 
is given by 

ry 1 4g 2 PW k 



H 3 n 



log i/. (5.6) 



Here, we have introduced the following expression for the logarithmic divergence, 

E E B2 = 4logiy - (5 - 7) 

u=-(j/-l)/2 R=\u/2\ 
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The wave function renormalization of the fermions can be computed in the same way [12]. 
The result is given by 

Next, we compute the renormalization of the coupling which can be read off from the trun- 
cated three point function of fermions and a scalar, 

We take the simplest choice of the external momenta: Jj = 0, Pi = 0, Sj = ij = (i = 1, 2, 3), 
mi — and K\ — K% — 1. In Fig. [71 all the diagrams contributing to the 1-loop corrections 
are listed. By following the same calculation that we described for the self-energy, we find 
the following values for the divergent part of each diagram, 

(i/-l)/2 n-l+u/2 „ 

(y-i) = 4,B OD (-ir + ' EE? 

^ u=-(u-X)/2 R=\u/2\ 

(v-l)/2 n-l+u/2 

(Y-2) = ^ ABCD (-l)^ E E 4 

^ u=-(i/-l)/2 fl=|«/2| 

(i/-l)/2 n-l+u/2 

ir-^Wir' E E 

" «=-(i/-l)/2 i?=|«/ 2 l 

(i/-l)/2 n-l+u/2 

(y-4) = A £ , BOfl( _ ir+ J E E 

p u=-(i/-l)/2 R=|m/2| 

(>-l)/2 n-l+u/2 , . 

(r-5)=4 K „i-ir> E E (-dp)- < 5 - 9 ) 

^ u=-(i/-l)/2 R=\u/2\ V 7 

Comparing the sum of the above diagrams with the tree level result, one obtains the renor- 
malization for the coupling at the 1-loop level, 

i 

This satisfies Z gpw = Z^Z£ and therefore the beta function is vanishing at the 1-loop level. 

6 Conclusion and discussion 

In this paper, we study a non-perturbative formulation of M = 4 SYM on R x S 3 proposed in 
[j~2] . After reviewing the formulation, we introduce Wilson loop operators in the formulation 
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(Y - 1) (Y - 2) (X - 3) 




(Y - 4) (Y - 5) 



Figure 7: The 1-loop corrections to the Yukawa coupling. The solid, dashed, curly and wavy 
lines represent the propagators of the 5*0(6) scalar, fermion, SO (3) scalar and the gauge 
field, respectively. 

which correspond to Wilson loops in M = 4 SYM studied in the context of the AdS/CFT 
correspondence. We calculate the VEV of the half-BPS Wilson loop by summing up all 
the planar ladder diagrams as done in the continuum theory in J39JS0] and reproduce the 
known result. We also calculate the one- loop beta function and verify that it vanishes in the 
continuum limit, which is consistent with restoration of the superconformal invariance. Our 
results serves as a check of the formulation. 

We should check that the planar non-ladder diagrams for the half-BPS Wilson loop are 
indeed canceled out at lower orders in the perturbative expansion, as done in [39], while 
the calculation in our case seems hard work. It is also desirable to derive the result using 
a method such as localization as done in the continuum theory in [15]. It is interesting to 
do the calculation of this paper in another non-perturbative formulation of M = 4 SYM on 
R x S 3 proposed in [36]. It is quite important to calculate the VEV of non-BPS Wilson 
loops such as a rectangular Wilson loop on R 4 in the strong coupling regime in the present 
formulation and reproduce the prediction from the gravity side [42] . This kind of calculation 
on the gauge theory side not done yet would serve as a highly nontrivial test of the AdS / CFT 
correspondence. 

We hope that our calculation in this paper will trigger development in the study of the 
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formulation. 



Note 



Preliminary result on the beta function in this paper was announced by G. I. and A. T. 
at the International Conference on Progress of String Theory and Quantum Fields Theory, 
Osaka City University 7-10 December 2007 [57l|58]. 
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A S 3 and the SU(2) group manifold 

In this appendix, we summarize some useful facts about S 3 and the SU (2) group manifold 
(see also p2}[15]). We regard S 3 as the SU{2) group manifold. We parameterize an element 
of SU (2) in terms of the Euler angles as 



where Jj are the generators of the SU(2) algebra, and < # < 7r, < ip < 2tt, < ip < 4tt. 
The periodicity with respect to these angle variables is expressed as 



The isometry of S 3 is SO (4) = SU(2) x SU(2), and these two SU (2)'s act on g from left 
and right, respectively. We construct the right-invariant 1-forms, 



9 = 



-ifJ3 -i9J2-ii>J3 



(Al) 



(9, <p, ip) ~ (9,<p + 2tt, V + 2tt) ~ (6, + 4tt). 



(A.2) 



dgg 



-i 



ifie l Ji 



(A.3) 
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where the radius of S 3 is 2/ fx. They are explicitly given by 



e 1 = — ( — sin ipd9 + sin 9 cos (pdip) , 
H 

e 2 = — (cos (pd9 + sin 9 sin cpdip) , 

e 3 = - (dip + cos 9dtp), (A.4) 



and satisfy the Maurer-Cartan equation 



2 

The metric is constructed from e* as 



de* - ^e ijk e j A e k = 0. (A.5) 



rfs 2 = eV = — (ci0 2 + sin 2 edip 2 + (di/> + cos 9d V ) 2 ) . (A.6) 

The Haar measure is defined through above metric as 

dQ 3 = - sin 9d9d<pdip, (A. 7) 

8 

which is left and right invariant. The Killing vectors dual to e l are given by 

d = --efd„ (A.8) 



where \i = 9,<p,ip and ef are inverse of el. It follows from ( fA~3l) and (lATSl) that 



£^ = -J^, (A.9) 

which indicates that Li are the generators of the left translation and satisfy 

[Ci,Cj] = ie ijk C k . (A. 10) 

The explicit form of the Killing vectors are 

C\ = —i ( — sin ipdg — cot 9 cos u?cL H ^cL ) , 

V smf / 

£ 2 = —i ( cos — cot # sin y?cL H %dii> ) , 

\ sin0 / 

C 3 = -id v . (A.ll) 
The spherical harmonics on S" 3 is defined through Wigner's D-function (see also [T3U50]): 



Y M (Q 3 ) = (-l) J - m V2JTTDj(-m,m) 
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= (-l) J - m V2J+l(J - m\g\Jm) 

= (-1) J ~"V2J+ le^-^iJ -m\e- iej2 \Jm), (A.12) 
where g is given in (1A.1I) . The complex conjugate is given by 

Y* mrh (ttz) = (-l) J - m v / 2TTTe- im0+ ^(J - m\e- l9j2 \ Jfh) 



= (-l) J - m V2J + le- im * +i ^(-l)- m -™(Jm\e- iej2 \ J - m) 

= {-l)"Yj_ m _^). (A.13) 

It follows from (IA.9[) that 

C±Y Jmrh (tt 3 ) = y> \J =F m)(J ± m + l)Y Jm±llfl , 

C 3 Y Jmfil = mY Jmih . (A. 14) 

By using the orthogonality relation 

J dg(Jm\g~ l \Jm){J'm'\g\J'm'} = ^j^^ JJ'Srnm'&mm, (A.15) 

where dg = dVL 3 /{2Ti 2 ) is the Haar measure, one can show the following equalities: 

^^ Y J imi m 1 (^) Y J2rn2m 2 {^3) = 8j lJ2 5 mi m 2 ^mim 2 (A.16) 



and 



(2J 2 + 1)(2J 3 + 1) 

2 J _|_ \ J2in 2 J 3 m 3 Ly J 2 m 2 Jzih 3 i / 



2 2 ^lWimj (fisjYjimzmQ (^3)^0377137713 (^3) 



I 1 \mi- m\r>J\— mi— mi f\ io\ 

— I X J U J 2 m 2 m 2 J 3 m 3 m3' I^.IOJ 

where is the Clebsch-Gordan coefficient. 

B Fuzzy spherical harmonics 

In this appendix, we summarize some useful properties of the fuzzy spherical harmonics 
(See H2HH], and also [3SIEIH5S])- 
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Given (2j + 1) x (2f + 1) rectangular complex matrices, we can generally express them 

as 



M = J2Mrr>\jr)(frl (B.l) 

r,r' 

where \jr) are the basis of the spin j representation of SU(2) algebra. {\jr){fr'\ | r = 
—j, ■ • • , j; r' — —j', ■ ■ ■ , j'} form a basis of (2j + 1) x (2f + 1) matrices. Let us consider 
linear maps, which map the set of (2j + 1) x (2/ + 1) complex matrices into itself, defined 
by 

U o \jr)(j'r'\ = Lf\jr)(j'r'\ - \jr){j'r'\Lf\ (B.2) 

where are the spin j representation matrices of the SU (2) generators. LjO satisfy the 
SU{2) algebra: [L^o, L^-o] =ie ijk L k o. 

By changing the basis, we can obtain the more appropriate basis for the action of L^o, 
which is called the fuzzy spherical harmonics: 

Yjmm = V^^-^'^rl'-rWij'r'l (B.3) 



where Cj*™ Jm is the Clebsch-Gordan coefficient and n is a positive constant, which is 
taken to be an integer as an ultraviolet cutoff in section 2. For a fixed J, the fuzzy spherical 
harmonics form the basis of the spin J irreducible representation of SU(2) under L^o, 

(Lio) 2 Y JmUj/) = J(J+ l)Vjm(jj')' 

£± ° ^JmOj') = a/(J =F m)(J ± m + ljyjmiiyjv), 

^3 ° ^Jm(ij') = mYjm^j,). (B.4) 

The hermitian conjugates of the fuzzy spherical harmonics are given by 

(*W)) f = (-l)" 1 -^^^)- ( B - 5 ) 

The orthonormality condition is 

-tr | ^JmOj')) = fijJ'Smm', (B.6) 
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where "tr" stands for the trace over (2f + 1) x (2f + 1) matrices. The trace of three fuzzy 
spherical harmonics is evaluated as 



AJimi(jj") _ ■ , i / , - \ v v • 

J2m 2 (jj')J3m 3 (j'j") ~ n } \ IJ l m l(M")j 1 J2m 2 (jj') 1 JsmsU'j") 



1 



(-l) Jl+i+i "y/<2J 2 + 1)(2 J 3 + 1)C£ 



■A ^2 
.hi»>-h.i»:i 1 »' j" j 



Oil 



(B.7) 



Cj 1 m 1 (j"j)J2m2(jj')J 3 m 3 j'j") — tr <j ^Jimi(j''j)^J 2 ra2(jj0^3 m 3(jY') 



^ ^ U J2m 2 (ji')^3m 3 (j'j")' 

where the last factor of (IB.7j) is the 6-j symbol. The 6-j symbol is related asymptotically to 
the Clebsch-Gordan coefficient: for R ^> 1, 

l _ 1 C_1 \-a+b+2c-d+e+2f+2R 

d + i? e + i? f + R)~ y/2R(2c+l) °— +/*-/+* 

We also introduce the vector fuzzy spherical harmonics Yj,^,^ and the spinor fuzzy 
spherical harmonics ^j m(j y) a , where p takes —1, 0, 1 and k takes —1 and 1. They are defined 
in terms of the scalar spherical harmonics as 



1 Jm(jj')i ~ 1 Z^ VmLj QplnQp(jf 



n,p 



Yj m (jj') a ~ ^2 C i;v7\cXup{jj') ' (B. 10) 

p 

where Q = J+8 p i, Q = J + 5 p -i and U = J + ^8 K i, U = J + V is an unitary matrix 

given by 

-10 1 

-i -i 
V2 



V = -^|-i (B.ll) 



The vector fuzzy spherical harmonics and the spinor fuzzy spherical harmonics satisfy 



ai) a pLi o +-S a j3 J KJ^y^ = k ( J + - J yj K m(ij , )Q . (B.12) 
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Their hermitian conjugates are given by 

(yp V — (_-\)m-(j-j')+i Y P 



Yjm(jj')a) — (-1)™ ^ j )+Ka+1 ^J- m ( 7 -' 7 ')-a, (B.13) 



Jm(jj')a J \ L ! 1 J-m(j'j)-ai 

and the orthonormality conditions are 



^tr j (YjmQj/yi^ Yj> m > (jj')ij ~ $J J'8mm'8pp> , 

^tr | (yj t m0y)a ) t ^ fc> J = 5/j'<W<W. (B.14) 

The trace of three fuzzy spherical harmonics, including the vector harmonics and/or the 
spinor harmonics, are evaluated as 



jyJm(j'j) ^ _1 ^ \ (y y Y Pl Y p2 

Jimi(j'j")pi Jimi{j" ])pi n \ \ J m ti'j) J Jimi(j'j")i J2rri2(j"j)i 

= \/3ii(2J + l)(2J t + l)(2.h + ip\ + 1)(2J 2 + 1)(2J 2 + 2p2 + 1) 

(-D^ +1+J ««'{! | , «/ *}. (B.15) 



X 



£jimi(jj')piJ2m2(j'j")p2J:im;i(j"j)p3 
1 



^6n(2J! + l)(2Ji + 2pf + 1)(2J 2 + 1)(2J 2 + 2p\ + 1)(2J 3 + 1)(2J 3 + 2p§ + 1) 
x( _ 1) -n±^^_o 1 -Q 2 -Q 3+ 2^W'(j|;l ^2 Qa ^ |0j Q 2 Qa\ (B16) 

IQs ft lj V m l ™2 m 3 / l J J J J 
"j) = ~ tr | (^jTmiO'j)a) ^J2 2 m 2 (j'j")a^JrnU"j) 



■pJimi(j' 'j)K\ 
J2m2(j'j")K2 Jm(j'' 



^2n{2lf 1 + 1)(2J + 1) 2 (2J 2 + 1)(2J 2 + 2) 
x | jJcS™ 1 Jm{? ? /'}■ (B.17) 



y J2m 2 (j'j")K 2 Jm(j"j)p - „ 11 I I 2 Jimi(j'j)a) ° c*P 1 J 2 m 2 (j'j")P 1 Jm(j"j)i 



= \/M5i + 1)(2J 2 + 1)(2 J 2 + 2)(2J + 1)(2 J + 2p 2 + 1) 

x (-l)f^wjS £ f}c£- Qm {? ^ «}, (B.18) 



where {• • • } with 9 slots is the 9-j symbol. 
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C Harmonic expansion of PWMM 



In this appendix, we make a harmonic expansion of the PWMM around the vacuum with 
(12.30!) an d (I2.33P for the perturbative analysis in section 4 and 5. 

The complete form of the action of the PWMM including fermionic part is given by 

Spw = -Jt- J dT Tr (^(^f + \( X i + leijk[Xj,Xk]) 2 + \d t <$> ab D t <$> ab + ^ab® AB 
- ±[X i7 <S> AB ][X t , <S> AB ] - ~[<S> AB , ® CD \[® AB , $ CD ] + + 3 ^ A ^ A 

where A, B are indices of 4 of SU(A). ip A and ip A are two-component spinors and $ab are 
5*0(6) scalars X m in (I2.56P rewritten in terms of SU(4) notation in the following way 

$ l4 = ^(X i+3 + iX l+6 ) (z = l,2,3), 

$ab = -$ba, ® AB = -® BA = & AB , ® AB = l -e ABCD <b CD . (C.2) 

For this expression, the harmonic expansion can be easily performed. 

By replacing Xi — > fiLi + Xi in (1C.1I) and adding the gauge fixing and the Fadeev-Popov 



terms 



Sgf+gh = -j— / drTr i i (d T A T + ifi[Li, Xi]) 2 - icd T D T c + /j l c[L i ,i/j l [L i , c] + i[c, X^] 1 
9pw J I 1 ) 



(C.3) 



then we obtain 

tgauge . ngauge . nmatter , nmatter 



q ngauge . ngauge . nmatter , nmatter (C A\ 

>~>PW+gf+gh — dpwjree D PW,int *~ ° PWJree D PW,int> V^'^) 



where 



S^free = ^ / *T Tr {^{d T X % f - £[L U A T f + \{d T A T f 

2 2 \ 

+!L(x i + ie ijk [L j ,X k ]) 2 - ^-{L^X,} 2 + icd 2 c-i^c[L i ,[L i ,c}}\ , (C.5) 

S P a w%t = 4- I dr Tr(i(d T X i )[A r ,X i ]+ f j,[A r ,X i ][L i ,A r ] + ^ r ,X 4 ] 2 
9pw J V ^ 

+ i mh {Xi + ie ilm [L h X m \)XiX k - -e ijk e ilm X j X k X l X m 
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~ i^[L u c\[c,X l } + d T c[A T ,c]j, (C.6) 

Sp^free = ^ J dr Tt (^d T <S> AB d T $ AB + AB - £[L U ® AB \[U, ^ B ] 

+^d T ^ A + ^\(^ A + a%^ A ])^j , (C.7) 

Spwtt = 4~ [ drTl Udr*AB)[Ar, ® AB ] - \[A T , $ AB ][A T , <S> AB ] - n[L h <S> AB ][X t , <S> AB ] 
9pw J \ A 

- l -[X^ AB \[X % M B ] - \[^ab^cd][^ AB ^ CD ] +t^ A [A T ^ A ] 
+^ A a i [X i ,i> A ] + i> A a 2 [$ AB , (4f] - ^ A ) T a^ AB ,^ B }) . (C.8) 

We perform a mode expansion for each (s, t) block for each field in terms of the fuzzy spherical 
harmonics defined in appendix A: 

js+jt J js+jt J 

J=\js-jt\m=-J J=\j s -j t \m=-J 

js+jt J js+jt J 

c (.,*) = ^ E c ^®W)> ^= E 

J=b's-jt| m=-J J=|j s -j t | m=-J 

</>*•<> =£ E E"<^'®^.,> 

K =±l(7 = | is _ it |m=-r/ 

-EE v-S"®«„o.«+ E E 

■Hjo-J'tl m=-J-| j=|j s _j t |_i m=-J 

v-r f =E E E^l®^,, 

«=±1 tr = |j s _j t | m=-I7 



j s +jt J +2 j s +jt-§ J 

k „/,(*> s )t ^ 

Jmijtjs) 1 



E E E E*£2~®V 

J=|j' s -J't| m=-J-i J=\js-3t\-h rn=-J 



1 js+jt Q 
X i ' = E E E ® ^Jm(i s j t )* 

' ,= - 1 Q=li s -it| m= -'3 

js+jt J+l js+jt J 

= £ £ ® ^™(jsjt)* + £ £ X ^™0 ® Y JmUsh)i 
J=\js-jt\ m=-J-l J=\j s -j t \ m=-J 
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+ ' if E 4£i ® ^i^. (C9) 

J"=lis— 1 m=-J 

Note that the modes in the right-hand sides are (2j s + 1) x (2j t + 1) matrices. The free part 
of (1C.4j) is expressed in terms of the modes as: 



ngauge . qmatter 
°PW,free <~ D PWJree 



+ \{-i) m - {j °- h) B^{-di+^j{j+ i)}Bf:i 

+ i(-l)^'-^K-MV(J + l)}c^i, 



(CIO) 



where 



ufj = fA(J+i), u + = rtJ + l), wj^^j + i). (C.ll) 

From (IC.lOp . the propagators can be read off as 

^(_l)m-( Js -i t ) +Mjj , 5m _ m , 5st ,^, 5 .^ ijfc27r5(p + (p = p' = 0) 

{B { si\v)ijB { s;£ ) { V ')ki) = ^(-l)^'- jt %j>6 m ^,6 gt ,6u>6 u 6 jk 2n6( P + p'^ ' 



n p^ + \i l J\J + l) 

(^(p^OOh) = ^(-l) m - (is - Jt) ^^5 m _ m ,^^^^ fe 27r5(p + p0 ^ 2+ * , 



{4>^AB Jm{p)io4>%B> J'm'{p)kl) = ~~^~7^^ABA'B'{ l) m {js Jt) 5jji5 m - m >5 s t>5ts>5ii5 jk 2lT5(p + p'^ 



2— jtl u/i y / U U "lib — ul v ol v uo v Lb v i r*, — ■ • ~ yj- i j- / 

p 2 + 

#PVKr r r s:A r r r r r, c/ A (*P + 



(^\p)i^iitAP%i) = —8jj>6 mn ,6 HH ,6i,6 a ,6 tv 6u6 s &c5(p - v V— J 2 - (C.12) 



" P 2 + c^ 



The gauge interaction terms in ( 1C.4I) are rewritten as 



ngauge 
J PW,int 
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n 
9pw 



dr tr 



lU J\m 1 {j s jt) J2m 2 (jtju)p2 J3m 3 (j u js)p3- U J 1 rn 1 \°T X J 2 m 2 p2 X J3m 3 p 3 ~ X J 2 m 2 p 2 ° tX J 3 m 3 p 3 ) 



A*( yJ^iJl + l)^Jimi(j s ji) J 2 m 2 (j t ju)0 J 3 m 3 (j u j s )p 3 \/Jl(Jl + 'j 2 m 2 (jtj u ) J3m 3 (j u j s )p 3 Jimi(j s j t )o) 

v R^'^ r(*> u ) ~.( u > s ) 
A ^ Jirri^ J 2 m 2 X J 3 m 3P3 

_l_ (_;n m -(js-ju)+i 

A V^Jirrnijsjt) J 2 m 2 (jtjv.)p 2 J-m(j u js)p ly J4m 4 (Jvjs) Jm(j s jt)p J3m 3 (j v .j v )p 3 1J J 1 -mi X J 2 m 2 p 2 x J 3 m 3 p 3 1J J 4 m 4 

— U Jim 1 (j s j t ) J 2 m 2 (j t j u )p 2 J-m(j u j s )p IJ J 3 m 3 (j u j v ) J 4 m 4 (j v j s )p 4 Jm(j s j u )p-D J im \ X J 2 m 2 p 2 J 3 m 3 X J 4 m4P4 > 

I • /Til \ p (s,*) (*,«) 

+ Wl\ J l + J-J^-JimiO'sjOpi J 2 m 2 (j t j u )P2 J3m 3 (j u js)p 3 X J imi pi X J 2 m 2 p 2 X J 3 m 3 p 3 

~ ' C J-m(j u js)p JxmxQsj^pi J 2 m 2 (j t ju)p2 C Jm(j s j u )p J 3 m 3 (j u j v )p 3 J4m4(j v js)p4 

(s,t) (t.u) (u,v) (v,s) 
Jim 1 p 1 ^ J 2 m 2 p2 J 3 m 3 p 3 J/ J 4 m 4 p4 

J 2 m,2 J 3 m 3 X J\m\p\ 
U J2m2{jsjt) J\m\{jt3u)p\ ■J3'm,3(jujs)0 c J 2 m 2 X J 1 m 1 p-i C J 3 m 3 ) 
~ ^Jim 1 (j s j t ) J 2 m 2 (jtj u ) J3rn 3 (j u j s )-Dj irni {0 T Cj 2rn2 Cj 3m3 + c J 2 rn 2 C> rC- J z rn 3 ) • (y>.Lo) 

The matter interaction terms in (IC.4I) are rewritten as 

5 matter 



PW,int 

n 



dr tr 



1 c ABCD/i , , , % R («,*) (Pt At,u) ,{u,s) _ ,(t,u) p. ,(u,s) > 

2 e u JimiOsit) J2m 2 Utju) J 3 m 3 (j u js) n J 1 m 1 l y T? ) AB,J 2 m2 ( PCD,J 3 m 3 V AB ,J 2 m 2 ar< Pc D,J 3 m 3 > 



1 1 ABCD/iJm(J s ju) n 

^ 2 Jimi(j«jt) J2m 2 (jtju) Jm (0s3u) J 3 m 3 (j u j v ) J 4 m4(j v j s ) 

v /dM tj{%u) ,(u,v) Ay,s) r(«,*) R («,f) \ 

X V- D Jimi- D J2m2^AB,J 3 m3^CD,J4m4 ^ Jimi ™AB, J2m2^J 3 m 3 ^CIV4m4 J 

m l(jsjt) J2iri2(jtjv.)0 J 3 m 3 (j u j a )p 3 
m 2(jtju) J 3 m 3 (j u j s )p 3 Jirni(j s jt)o) ( t > AB,Jimi ( l ) CD,J2rn2 X J 3 m 3 p 3 



^ ABCD(_ 1)m -(j s -j u )+l 
2 v 7 



l^J4m 4 (>j s ) Jm(j s ju)p J 3 m 3 (jujv)p 3 Ly J2m 2 Utju) J-m(j u j s )p Jim 1 (j 3 jt)pi x J 1 m 1 p 1 l f J AB,J2m2 J 3 m 3 p 3 rCD,J4rri4 

_ -fy -I s ) ( s >*) aS U,V ^ ( V,S ^ \ 

u J4rri4(j u j v ) J 3 m 3 (j v j s )p 3 Jm(j s j u )p u J2m2(jtju) J-m(j u j s )p Jim 1 (j s j t )pi x J l m 1 p 1 'rAB,J2m2 ( PCD,J4m4 X J 3 m 3 p 3 ) 

__ 1 ABEF CDGHnJm{j s j u ) A 

(jsjt) J2m 2 (jtju) Jm tisju) J 3 m 3 (j u jv) J4m 4 (j v j s ) 

X f^ (s,t) r/) (n, '' ) rii^'^ - rti (s,t) r/)^'^ ^ 

A V^AB, Jimi ^CD,J 2 rn 2 ^EF,J 3 m 3 ^GH, J 4 m 4 ^AB,J imi ^EF, J 2 m 2 ^CD,J 3 m 3 y j GH,J 4 m 4 > 

+ i f (_ 1 )r^-(j.-i„)+^^a-^(jti»)« w,(;.t)t ^(«,t) 

\ Ji-mi(j t j s )«i J 3 m 3 (j s j u ) ^A,Jimi/ti J 3 m 3 ^J 2 m 2 K 2 

_-t r Jim 1 (j s j t )K 1 il>( s,t ^ / „(u,t) \ 

■hm2{j B 3 v )K2 J 3 m 3 (j u jt)^A,J 1 m 1 K 1 Y J 2 m 2 n 2 J 3 m 3 J 

_ ( (_\\m 3 -{js-ju)+ !: ^^QJ2-rn 2 {jtju)K 2 ,{s,t)] (s,u) ,A(u,t) 

\} 1 y Ji-m.i(jtj s )fti J 3 m3(j a ju)p3 rA,JxmiKi J 3 m 3 p 3 ^ J 2 m 2 K 2 
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, ( jJim 1 (j a j t )K 1 / (s,*)t jMa,u) (u,i) \ 

~ ry J2m2(j s ju)K2 J 3 m 3 (j u j t )p 3 i f J A,J 1 m 1 K 1 r J 2 ma«2 J 3 m 3 p 3 J 

% _ABCD ( ( _,s mi -(j s -j t )-?±fJ2m2(Jtju)K2 ,{s,t)] ,(s,u) , (*,«)t 

2 V ' ^Ji-m 1 {j t j a )K 1 J 3m3 {j s j u )^A,J imi K 1 ^CD,J 3m3 ^B,J 2 m 2 K 2 

,(_-,\m 2 -(ju-js)-^ 2 2 -fJinii(jsjt)Ki 7/,( s '*)t ,/,(«. s)t \ 

" 1 ; J 2 -m 2 (j 3 j u )K 2 J3m3(j I1 j t )^A,J 1 m 1 K 1 ^B,J2m 2 K2^CD,J3m 3 y 

j ( (_l)-m2-(ju-js)-^-f J 2-m 2 (jsju)K2 ,A,(s,t) ,(t,u) j^B(u,s) 

\y ' Jimi(j s jt)Ki J 3 m 3 (jtj u ) ^ J\m\K\r AB,J 3 m 3 ^J 2 m 2 K 2 

_(_|'_1 N | m i-(is-ii)-^ fJi-miUtjs)^ ,i,Ms,t) jB{t,u) ,{u,s) \ 

' J 2 m 2 (jij u ) K2 J 3 m30 u i s )^-/im 1 Ki ( / y J 2 m 2 K 2 ^A_B,J3m 3 y 



(C.14) 



D Generalization of the Feynman formula 



We first show that for pi , • • ■ , p k > 

dri ■ ■ ■ / dr k 8{l - n r fc) T f 1_1 ■ ■ • r > P ^ 



! = r( Pl )---r( Pfc ) 



(D.l) 



'o J o 

The left-hand side of (ID.lj) equals 

'1 /•!— Tl yl— TJ T fc _ 2 



ii Jo Jo 



Tfc-i)'*- 1 . (D.2) 



By putting Tfc_i = (1 — T\ — • • • — Tfc_ 2 )<7fc_i, we calculate this as 

'1 fl—Ti n T fc _ 3 



P(p k - U p k ) / dri / dr 2 --- 

/o Jo Jo 



«^fc-2^1 T k-2 l 1 ~~ T l 



By repeating this procedure, we see that the left-hand side of fID.ll) equals 



/3(p k -i,Pk)/3(Pk-2,Pk-i +Pk) • ■ -P(pi,P2 H I-PjO- 



By using 



/3(a,6) 



r(a)r(6) 



r(a + b) ' 

we see that this is indeed equal to the right-hand side of (ID. II) . 
Here we would like to show that 

1 r\ m 

-e A +a B 

m 

9=0 



m! dg r 



Tk-2 



jPfc-i+Pfe-1 



(D.3) 



(D.4) 



(D.5) 
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= / cbn--- I da m+1 S(l-a! a m+1 )e aiA Be a2A B ■ ■ ■ Be amA Be a ™ +lA 

Jo Jo 

pi /"l-ai pl—ai a-m-i 

= da x da 2 --- da m e aiA Be a2A B ■ ■ ■ Be amA Be (1 - ai —~ am)A . (D.6) 

Jo Jo Jo 

By expanding the exponential and using (ID. II) . we calculate the second line of (ID. 61) as 



oo If 1 f 1 

y2 — j da!--- da m+1 5(l-ai a m+ i)a'/ • • • a|^J 

, , n n- ■ ■ ■ Wi ! Jo Jo 

x A h BA h B ■ ■ ■ BA lm BA lm+1 

= Y — nh + 1) • • • r (Wi + 1) . . . BA lm BA lm+1 

, / 1 !---/ m+1 !r(/ 1 + ••• + /„+!+ m + l) 

X ^ 

V -rA ll BA h B ■ ■ ■ BA lm BA lm+1 

, n (Zi + --- + Z m+1 +m ! 



X 

= ^ A h BA h B ■ ■ ■ BA lm BA lm+1 . (D.7) 

n=m iiH h£ m +i=n-»Ti, Zi,--- ,Z m+ i>0 

This agrees with the left-hand side of ( ID. 61) . 
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